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1 Introduction 



Recently there have been attempts to derive quantum dynamics (canonical 
commutation relations and Heisenberg equations of motion) as an emergent 
phenomenon. Interesting results in this problem are obtained by S. Adler and 
collaborators in a series of papers Pfl El EH E] where the quantum dynamics 
is derived, under certain approximations and assumptions, from the classical 
statistical mechanics of matrix models. 

Here we want to present a result which might be useful in clarifying 
certain aspects of this emergent phenomenon. We consider the statistical 
mechanics of a + 1 dimensional Hermitean matrix model in the presence 
of a random external source with a Gaussian distribution and show that in 
the high temperature limit the conjugate variables (conjugate with respect 
to the Poisson brackets) obey the canonical commutation relations. Since 
this result is obtained in the high temperature limit it is independent of the 
specific form of the Hamiltonian of the system. Thus the only assumption we 
make about the structure of the matrix models is the randomness of the ex- 
ternal source. We derive our result by using Witten's nonabelian localization 
principle 2 jHJ El 13 El • We show that the term in the partition function which 
survives the high temperature limit is a closed equivariant form. Then we 
use the nonabelian localization principle to show that the partition function 
is localized on a certain subset of the phase space. The subset in question 
is characterized by a 2 dimensional Yang-Mills algebra. This algebra is an 
example of a cubic Artin-Schelter algebra and is the universal enveloping 
algebra of the Heisenberg algebra. 

In Sec. 2 we construct the phase spaces of matrix models that we will 
consider. We show that the unitary group U(N) has a Hamiltonian action on 
the phase space and derive the corresponding moment map which is conserved 
whenever the Hamiltonian is U(N) invariant. 

In Sec. 3 we work out the statistical mechanics of matrix models defined 

2 The adjective nonabelain refers to the noncommutativity of a symmetry group. The 
abelian localization was derived by Duistermaat and Heckman 
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in Sec. 2. Since the moment map is a conserved quantity for U(N) invariant 
systems we expect the Boltzmann factor to depend not only on the the 
Hamiltonian but also on the moment map. This dependence on the moment 
map is responsible for a nontrivial high temperature limit of the partition 
function. We assume that the external source coupling to the moment map 
is itself random with a Gaussian distribution (white noise). This assumption 
allows us to apply the nonabelian localization to the partition function in 
the high temperature limit. 

In Sec. 4 we review the basics of equivariant cohomology and nonabelian 
localization principle. 

Finally in Sec. 5 we apply the nonabelian localization principle to the 
partition function derived in Sec. 3 and show that the latter is localized 
on a subset of the phase space characterized by the canonical commutation 
relations. 

2 Matrix Models 

A matrix model is a dynamical system whose configuration space is a subset 
of the space of N x N real or complex matrices. In most applications the 
configuration space is either a vector space (for example hermitean matrix 
models) or a group manifold (for example unitary matrix models). In this 
article we will work with hermitean matrix models. Thus the configuration 
space C will be the space of all N x N hermitean matrices. The phase space 
Jv[ = T*C of the system will be identified with C 2 . The coordinates in M. 
will be denoted by (X, P). We will also use the notation X 1 = ReX and 
X 2 = ImX, and similarly for P. 

A symplectic form on M. is defined as follows 

2 

u = dX rij A dP rij = Re Tr (dX A dP f ). (1) 

r=l 

The action of the unitary group U (N) given by 

X -> UXU' 1 (2) 
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P -> UPU- 1 (3) 

is symplectic. 

The action of the fundamental vector field corresponding to an element 
— iA (A 1 " = A) of u(N) is given by 

S A X = i[A,X] (4) 

5 A P = i[A,P}. (5) 

In terms of real and imaginary components this is equivalent to 

5 A X 1 = -[A,X] 2 (6) 

5 A X 2 = [A, X] l . (7) 

and similarly for P. Thus the fundamental vector field itself is given by 



Proposition 2.1 The action ofU(N) on M. is Hamiltonian. 

Proof: We must show that (i) the fundamental vector fields are Hamilto- 
nian and (ii) that the corresponding moment maps form a Lie algebra (with 
respect to the Poisson bracket induced by uo) homomorphic to Q 

(i) The fundamental vector fields are Hamiltonian i.e. there exist a func- 
tion ji A on M. such that 

i u u = e rs [A, X] rij dP sij - e rs [A, P\ Hj dX sij = -dfi A - (9) 

Let 

H A = Tr(-iA[X,P}). (10) 

then 

Tr(-iA[X,P}) = ImTr (A[X,P}) (11) 
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= - ImTr ([A,P]X) (12) 

= ImTr ([A,X]P) (13) 

but 

ImTr([A,P]X) = Tr ([A, P] 1 X 2 + [A, P] 2 X±) (14) 

= - e rs [A, P] rij X sij (15) 

thus 

°" A =-e r .[A,P] rii . (16) 



Similarly, 

-^- = trs[A,X] rij . (17) 

So we conclude 

- d/i A = i u u. (18) 

(ii) Let -%T a be a basis for u{N) with Tj = T a , [T a , T b ] = if ahc T c and 
Tr T a T 5 = \5 ab . 
Define 

£a = ^T a = e rs ([T a ,X] r ..^ + [T a ,PU-|-] (19) 
ix a = ^ Ta =Tr(T a [X,P]). (20) 
We want to show that the action of u(N) is equivariant i.e. 

{fJ-a, fJ-b} = ifabc^c- (21) 

Here { , } is the Poisson bracket corresponding to the symplectic form uj. 
We have 

{/ia,/ifc} = w(£o,&) (22) 
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= (<€.<")(&) (23) 

= Tr([T a ,P][X,T b ] + [T a ,X][P,T b ]) (24) 

= Tr[T a ,T b ][X,P] (25) 

= ifabcHc- (26) 



Hence we conclude that the action of u(N) on M is Hamiltonian. 

Now that we have the phase space, our next task is to choose a Hamilto- 
nian. Matrix models whose Hamiltonians are given by traces of polynomials 
in X and P are the most common ones in physics and geometry. Such trace 
Hamiltonians are invariant under the action of U(N) and consequently the 
moment maps \x a are conserved quantities. In what follows we will assume 
that the dynamics is U(N) invariant; thus the Hamiltonian will be taken as 
a trace Hamiltonian. 

Finally let us note that the Liouville measure on M. is given by 

— = du(X)du(P) (27) 
nl 

where 

du{X) = HdX u J] J[dX rij (28) 

i r=l,2i<j 

and similarly for du(P). 

3 Statistical Mechanics of Matrix Models 

For a system in equilibrium the distribution function p of the statistical 
ensemble should be a conserved quantity. Therefore one expects p to be a 
function of the conserved quantities of the system. So in the case of matrix 
models we have 

p = p(p a ,H). (29) 
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The explicit form of the equilibrium distribution can be derived by maximiz- 
ing the entropy 

S = - f^-r pin p (30) 
J n\ 

while keeping the mean values of 1 (normalization), p a and H fixed. Thus 
effectively one maximizes the functional 

_ p6 _j _ p4 ,y»-J- pP H (3i) 

where 9, <p a and (3 are the Lagrange multipliers. The result unit normalization 
is 

p = L e -n*-f* (32) 

Zj 

where Z is the partition function 

Z = [ ^ e -<t>>«-PH. (33) 



The same result can also be derived by counting microstates. This alternative 
derivation is given by Adler and Horwitz in 

Using the explicit form of the moment maps and defining the hermitean 
matrix = <j) a T a we have 

Z = J e "-^[x,P]-pH_ (34) 

As we will see in the next section the combination 

to = uj - (j) a p a =u- Tr<f>[X, P] (35) 

appearing in the partition function has an important geometric meaning and 
is called the equivariant symplectic form. 

From a field theoretic point of view the Lagrange multiplier cf> plays the 
role of an external source coupled to a conserved charge of the system. In 
our model we will assume that this external source is itself random with a 
Gaussian distribution 3 . Thus the modified partition function of the system 



5 A more familiar example of a situation where the source term is random would be a 
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with random external source is 

(36) 
(37) 

4 Equivariant Cohomology and Nonabelain 
localization 

In this section we review the basics of Witten's nonabelain localization prin- 
ciple. Let (A4,cu) be a symplectic manifold and G a compact Lie group with 
a semi-simple Lie algebra Q. Suppose G acts on M. with a Hamiltonian ac- 
tion i.e. assume that the action is symplectic, the fundamental vector fields 
are Hamiltonian and the moment maps form a Lie algebra (with respect to 
the Poisson bracket) homomorphic to Q. The fundamental vector field corre- 
sponding to an element £ G Q will be denoted by V(£). If {e a } is a basis for 
Q then the dual basis will be denoted by {e a } and as a short hand notation 
V a = V(e a ). Similarly the Lie derivative along the direction of V a and the 
contraction with V a will be denoted by £ a and i a respectively. 

Consider the tensor product Q(Ai) ® P(G*) where Q(A4) is the exterior 
algebra on Ai and P(Q*) is the polynomial algebra on Q. A generic element 
in Q(M) <g> P(Q*) is of the form 

E ( 38 ) 

k 

where a k is a k-form on M. and Pk(4>) is a polynomial in Q. The elements of 
magnetic system coupled to a random external magnetic field. 



e ZJ-f3H-eTr4> 2 



where 

is the Gaussian distribution of the source (p. Note that Tr<fi 2 
the Cartan-Killing form on u(N). 



7 



Q(M) <S> P{G*) can be multiplied as 

E «*W)E PQitt) = EE(« fe A^)P fc (0)QK0) (39) 



Since any polynomial is a linear combination of monomials a generic ele- 
ment in Vt(M)®P{Q*) can be written as a linear combination of the elements 
of the form 

a p <p ai ...0 a \ (40) 
This also defines a grading in Q(M) <E> -P(£/*) 

fi(.M)®P(S*) = (41) 

A fc A z C A k+l (42) 

A fe = span {a p (j) ai . . . ^ : p + 2r = k} (43) 

With this grading the multiplication is supercommutative. 
There is an action of Q on Q(A4) <g> P(Q*) given by 

£t <g> 1 + 1 <g> Q (44) 

where £g is the Lie derivative along the fundamental vector field V(£) and 
is the coadjoint action of Q on P(Q*) given by 



r 



C t {<t> ai . . . 4> ar ) = E ai • • • a "" 1 (-rcX0 6 )0 a " +1 "'* Or - (45) 



71=1 



Here c c ab are the structure constants of Q. Also note that C l j ) ((f> ai . . . <p ar ) = 0. 
The Cartan differential on Q(M) <S> P{Q*) is defined as 

D = d <g> 1 - i a <g> a . (46) 

In accordance with the isomorphism P(Q*) = <S'(^) we assume that under a 
change of basis of £ the variables (f> a transform as the elements of the dual 
basis. Then it is easy to see that the combination i a ®(p a is basis independent. 
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It is also straightforward to show that D is a superderivation of degree 1 

DA k c A k+l (47) 

D(ab) = {Da)b+{-l) k a(Db) (48) 

for a G A k . 
Proposition 4.1 

D 2 = —£ a <S> <p a (49) 

Proof: 

D 2 = (d(g)l-i a (g)(f) a )(d(g)l-i b (g)(j) b ) (50) 

= d 2 ®l- di a ®(j) a -i a d®4> a + i a i b ®4> a (j) b (51) 

= ~£a®r (52) 

The G-invariant part (Q(M) P{Q*)) G of Q(M) <g> P((?*) is called the 
Cartan complex and is denoted by Qg{M)- 111 other words Qg{M) is the set 
of all elements of Q(M) <E> P{Q*) which get annihilated by <g> 1 + 1 <g) Q 
for all £ G 0. 

Proposition 4.2 

(A <g> 1 + 1 <g> C a )£> = D(£ a <g> 1 + 1 <g> C„) (53) 

Proof: 

[£ a ®l + l®C a ,P] = -[£ a ,i 6 ]®0 6 -i 6 (g)(C a 6 -0 6 C a ) (54) 

= ® 6 - ® (-40 c ) (55) 

= -c^(i c ®0 6 + i c ®(c^ 6 ) = O (56) 

As an immediate consequence of this result we see that the Cartan com- 
plex inherits the grading in Q(M) ® P(Q*) 

MM)®p(g*)f = ®A G (57) 
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where is the G-invariant part of A/.. 
Proposition 4.3 On Q G (M) 

D 2 = 0. (58) 

Proof: Let b G &,g(M) then 

(£ a ®l + l®C a )6 = => (l®0 a )(£ a ® l + l®C a )6 = 0^ (59) 

(£ a ®0 a + l®0 a C a )6 = O^ (60) 

(-D 2 + 1 ® a C a )b = (61) 

but as we noted earlier (1 ® (f) a C a )P((f)) = 0. So G-invariance of implies 
D 2 = on fi G (M). 

This result allows us to construct the cohomology of D on the Cartan 
complex. The resulting cohomology is called the equivariant cohomology of 
M. with respect to the group G and is denoted by Hq(A4). 

Now we have to show that all this construction is not vacuous by giving 
an example of an equivariant form. The form of degree 2 defined by 

ZJ = lu - (j) a /j a . (62) 

is called the equivariant symplectic form. Here one can easily show that the 
term <fi a fi a is basis independent. It is again easy to show [T2| that ZJ is G 
invariant and Du = 0. This equivariant symplectic form is encountered quite 
frequently in applications of equivariant localization principle to problems in 
physics and geometry (see for instance El EH!)- 
Equivariant forms can be integrated as 

/ J2 * k Pk(<P)= I d<l> l ...d<t> N e-^P k {<ft I a k (63) 
J ^ Jq Jm 

where ( , ) is the Cartan-Killing form on Q and is assumed to be nondegen- 
erate i.e. Q is assumed to be semi-simple. The exponential factor e~ € ^'^ 
is introduced to the definition as a convergence factor for the integral over 
the Lie algebra. Note that e(0, 0) is an equivariant form of degree 4. The 
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integral of a k over M. vanishes unless k = dim M.. 



Proposition 4.4 If M. is without boundary then the integral of any D-exact 
form vanishes. 



Proof: Note that 

D J2 = E {du k )P k {4>) - {i a a k )4> a P k {4>)- 



(64) 



Then the integral of the first term over M. vanishes because of Stokes' the- 
orem, while the integral of the second term vanishes since exterior degree of 
i a a k is always less than dim M.. 

Proposition 4.5 Let a be a closed equivariant form and X be an equivariant 
1-form. Then 



I a = I a 



.itDX 



(65) 



Proof: The proof is based on the observation that the difference of the two 
integrands is D-exact 



AtD\\ 



= (1 


- e ltDX )a 


(66) 


= -(UDX+ UitfDXD\+ ...)a 


(67) 


= D 


r l 

-(itX + -(it) 2 + XDX+ . . > 


(68) 


+(itX + -(it) 2 + XDX + . . .)Da 
2 


(69) 


= D 


-(itX + hit) 2 XDX + . . > . 


(70) 



Now the desired result follows as a simple corollary to the previous proposi- 
tion. 

The localization principle follows from the evaluation of the integral 



J a e UDX 



a e itd\ e -it^\{Va) 



(71) 
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in the asymptotic limit t — > oo. In this limit the integral can be calculated 
by the method of stationary phase which implies the localization of the 
integral on the critical points of the function 

a A(K). (72) 
There are two sets of equations for the critical points: 

A(K) = (73) 
(P a dX(V a ) = 0. (74) 

In order to ensure the compactness of the critical point set we will assume 
that the last equation has the unique solution (j) a = 0. 

In general there is no restriction in the choice of the equivariant 1-form A. 
We can use this freedom to our advantage and choose a convenient A whose 
critical points are (relatively) easy to find. In most physical and geometric 
applications such a convenient choice is given by 

A = Jdl. (75) 

Here J is a G-invariant almost complex structure compatible with the sym- 
plectic form u: £ a J = 0, u(JX, JY) = w(X,Y) and g(X,Y) := u(X, JY) 
is a positive definite metric on M; and 

I = E& ( 76 ) 

a 

Proposition 4.6 A is G invariant. 
Proof: 

£ a JdI = (£ a J)dI + 2J£ a (fi b dfi b ). (77) 

The first term on the right hand side vanishes since J is assumed to be a 
G-invariant. For the second term on the right hand side we have 

2J£ a (ij, b d^ b ) = 2J((£ a fi b )d{i b + /2 b d(£ a ij, b ) (78) 

= 2J({fi a ,fi b }dfi b + [i b d {n a , fi b } (79) 
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2Jd({fx a ,n b } dji b ) 



(80) 



UdiifabclXcUb) = 



(81) 



since the structure constants ifabc of a compact, semi-simple group are com- 
pletely antisymmetric. 

The following result is very important in the practical applications of 
nonabelian localization. Its proof is based on the compatibility of J with u 
and is given in 

Proposition 4.7 The critical points of the function (p a X(V a ) coincide with 
the critical points of the function I. 

Thus there are two types of critical points: 

• Ordinary critical points p: /U a (p) = for all a 

• Higher critical points p: dl p = but not all /i a (p)'s vanish. 

The localization principle can now be applied to the integral of e u . 



and implies that the integral is localized on the critical point set of /. 

5 Localization for Matrix Models 

Now consider the partition function of the Hermitean matrix model in the 
high temparature limit j3 — > 



is nothing but the Cartan-Killing form of u(N). This integral is clearly of 
the form discussed in the previous section. Therefore we conclude that the 




(82) 




(83) 



But as we noted earlier the term 



Trcf) 2 



(84) 
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integral is localized on the critical point set of 



Y^Tr (T a [X,P])Tr (T a [X,P]) 



(85) 



a 



a 



(86) 



The critical points are given by the equation 



dI = Tr[ P, [X, P] ]dX -Tr[X, [X, P] ]dP = 



(87) 



which implies 



[ P, [X, P] ] 







(88) 



[ x, [X, P] ] 



0. 



(89) 



Let At = X, A 2 = P and [A^ ® A v ] = A^ ® A v - A v ® A^. 

Now the tensor algebra generated by span {Ax, A 2 } modulo the two sided 
ideal generated by 



is an example of a Yang-Mills algebra in 2 (euclidean) dimensions. This 
2 dimensional Yang-Mills algebra is also an Artin-Schecter algebra and is 
known to be isomorphic to the universal enveloping algebra of the Heisenberg 
algebra |TH ITol ITol 117] . In particular this means that the only solution for 
dl = is the canonical commutation relation 



where a is a real number and 1 is the unit matrix. 

Thus we conclude that at high temperatures the coordinates of the phase 
space of a Hermitean matrix model obey the canonical commutation rela- 
tions. 

The ordinary critical points correspond to the commutative case a = 0. 
In this case localization is onto the Cartan subalgebra of u(N). In fact at 
finite iV the ordinary critical points are the only critical points since there are 




(90) 



[X,P] 



ial 



(91) 
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no finite dimensional representations of the Heisenberg algebra. The higher 
critical points (a 7^ 0) exist only in the large N limit. 
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